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Summary. We use Euler's difference lemma to prove that, for > and < A < 1, the 
function P n defined on the non-negative integers by 

defines a probability distribution, known as the Generalized Poisson Distribution. 
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1. Introduction 

The Generalized Pois son Dist r ibutio n (GPD), introduced in IConsul and Jain] l)l973|) . and 
studied extensively bv IConsul l)l989(l is defined on the non-negative integers, for < A < 1 
and > 0, by 

n\ 

Applications of the GPD can be found in settings where one seeks to describe the 
distribution of an event that occurs rarely in a short period, but where we observe the 
frequency of its occurrence in longer periods of time. It extends the Poisson distribution by 
its ability to describe situations where the probability of occurrence of a single event does 
not remain constant (as in a Poisson proce ss) , but is affected by previous occurrences. The 
distribution has been found l|Consu]Lll989t pp. 117-129) to accurately describe phenomena 
as diverse as the observed number of industrial accidents and injuries, where a learning 
effect may be present, the spatial distribution of insects, where initial occupation of a spot 
by a member of the species has an influence on the attractiveness of the spot to other 
members of the species, and the number of units of different commodities purchased by 
consumers, where current sales have an impact on the level of subsequent sales through 
repeat purchases. 

It was shown bv lConsul and Jairl l|l973j) that JIJ is a probability distrib ution, as it has 
the property YlnLo Pn(0 > — 1 , by using an identity that can be found in IJensenl (l902, 
eq. 6). However, as remarked bv lConsul l)l989t p. 12), "It is very diffic ult to prove by direct 
summation that the sum of all the probabilities is unity" . Recently, iLerner et al J Jl997l) 
gave a more direct proof using analytic functions. We give a shorter and more elegant proof 
based upon an application of Euler's classic difference lemma. 

2. Derivation 

We shall first prove the identity 

£ ^±^e-^« = JL-, for - Ao < A < 1, (2) 
n\ 1 — A 

71=0 
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where A = 0.2784645428 ... is the solution to Ae A = e _1 . Let S(0, A) denote the sum in 10. 
Expanding the exponential, grouping terms and changing the order of summation gives 

S(M) . f^f^ff^ . £ B-D-f^ (3) 
= E^E(-Wf> + ^. ( 4) 



fc! 

fc=0 n=0 



Now use Euler's difference formula, which states that the fcth difference of any pth power 
is for p < k and fc! times the leading coefficient of the pth power for p = fc. The most 
convenient representation for our purpose, as can be found in iGouldl l)l978l eq. 5.12), is 

k \ , a , o^p _ ( 0. 0< P <k 



n=0 



Applying this to the inner summation in (@J gives S(8, A) = If^ ^' = Sfclo ^ fe = r=A > 

and establishes the identity in J2J. 

We now show that for |A| < Ao the interchanging of the summation signs in the proof 
of (J2J is allowed, as we have tacitly assumed, by establishing absolute convergence. Taking 
absolute values of the summands one sees that the inner summation in reduces to el e+A " . 
Now apply Cauchy's root test and Stirling's approximation n\ ~ n n e~ n \/2im to give 



l imsup ? \ e + M n e \6+\n\ = |A|e 1+ l A l < 1, 
rwoo V n\ 

or the desired |A| < Ao as a criterion for absolute convergence. Another application of the 
root test shows that the left-hand side of J2J converges for |Ae~ A | < e" 1 , and as all the 
summands are positive, this convergence is uniform. Thus we can extend the range by 
analyticity to — Ao < A < 1, and this completes the proof. 
Using formula J2J, we have 

^— ' ^-^ nl (n—ly. 

n=0 n=0 n=l v ' 



-e 



_+\ + n\) n n _ e _ x _ nX 

n=0 



= S(6,\) - AS(6» + A,A) = 1, 
and proves that Q indeed defines a probability distribution 
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